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ABSTRACT. The Ricci flow is a parabolic evolution equation in the space of Rie- 
mannian metrics of a smooth manifold. To some extent, Einstein equations give rise 
to a similar hyperbolic evolution. The present text is an introductory exposition to 
Bianchi-Ricci and Bianchi-Einstein flows, that is, the restricted finitely dimensional 
dynamical systems, obtained by considering homogeneous metrics. 

1. Introduction 
These notes are variations around the homogeneous space 

X n = Symt = GL(n,R)/0(n) 

that is, the space ofnxn positive definite symmetric matrices, and sometimes, its 
subspace of those matrices with determinant 1, 

Y n = SSym+ = SL(n,M.)/SO{n) 

Besides their striking beauty, these spaces modelize many structures and support 
fascinating geometry and dynamics. It is surely very interesting to investigate in- 
terplays between these aspects. We will not do it systematically here, but rather 
briefly note some of themj. Our modest remark here is that "Bianchi" Ricci flows 
and Bianchi cosmologies are better seen as natural dynamical systems (i.e. differen- 
tial equations) on X n , respectively of first and second orders. To be more precise, 
Bianchi spaces are special homogeneous Riemannian spaces, those given by left in- 
variant metrics on Lie groups. For a given Lie group G, the space of such metrics 
is identified with X n , n = dimC7. The Ricci flow acting of the space of Riemann- 
ian metrics on a manifold, becomes here a (gradient-like) flow on X n . Similarly, 
the Cauchy problem for the (vaccum) Einstein equations becomes here a dynamical 
system on TX n . The group G acts preserving all these dynamical systems. 



We would like to emphasis on that this is a preliminary short non-finished work. This explains 
in particular why many proofs are left as exercises. Also, this text may be considered as expository, 
although the method here does not follow any existing approach. 
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References. There are now abundant references on the Ricci flow techniques, since 
their use by Perelman in his program on the geometrization conjecture on 3-manifolds. 
We may quote as an example iflOl as an interesting recent reference. The "toy" ho- 
mogeneous case (with which we are dealing here) was in particular investigated in 

nasi. 

The history of cosmology from a relativistic point of view, i.e. applying Einstein 
equations, is very old, and was in fact usually considered within a homogeneous 
framework, even an isotropic one, as in the standard big -bang model lfl2l l24l l25ll . 
Even, concrete interplay between cosmology and the mathematical theory of dynam- 
ical systems, were involved in the literature, but this seems to be not well known by 
"mathematicians" (which gives motivation for our text here). As recent references, 
we may quote |[2l l2Tll20ll . 

For the general material on homogeneous spaces, we quote ll6ll9l [T7l[T8l . 

Finally, in these proceedings, recommended references would include Ifl5ll22l . 

2. A MULTIFACETED SPACE 

The general linear group GL(n, R) acts transitively on the space of positive scalar 
products on R n , the stabilizer of the canonical scalar product being the orthogonal 
group 0(n). This space is therefore identified to the homogeneous space GL(n, R) /0(n). 
It will be sometimes more convenient to deal with a reduced variant: 

Y n = SL(n,R)/SO{n) 

which then represents conformal scalar products, or equivalently, scalar products with 
unit volume (i.e. their unit ball has volume 1 with respect to the canonical volume). 

For a more intrinsic treatment, we start with a real vector space E, and consider 
Sym(E) the space of its quadratic forms. Inside it, we have the open subspace of all 
non-degenerate ones Sym*(E), and Sym + (E) (or X(E)) the open cone of positive 
definite ones. We also get spaces of conformal structures by taking quotient by R 
acting by homothety; in particular S Sym + (E) (or Y(E)) will denote the space of 
conformal positive definite structures. In the case of E = R n , we use the notations: 

Sym n , Sym* n , Syrn+(= X n ), SSym+(= Y n ) 

They are identified to subspaces of symmetric matrices {A = A* G M nn }. The 
last three spaces correspond respectively to: det A / 0, A has positive eigenvalues, 
and positive eigenvalues with det A = 1. 

Exercise 2.1. Show that the GL(n, M.)-action on Sym n is given by g.A = g* Ag 
(where g* is the transpose of g). 

2.0.1. A metric on the space of metrics. A Euclidean structure q on a vector space 
E induces similar ones on associated spaces, in particular on the dual E* and on 
E* ®E* . If (ej) is a g-orthonormal basis, then its dual basis (e*) is also orthonormal, 
and also is the basis (e* (g> e*). 



HOMOGENEOUS SPACES, DYNAMICS, COSMOLOGY.. 



3 



Now, since Sym+ is open in Sym n , its tangent space at any point q is naturally 
identified to Sym n , which is thus endowed with the scalar product (,) q . Therefore, 
Sym+ becomes (tautologically) a Riemannian space. 

Exercise 2.2. Show that 

(1) {p,p) q = tr^pq^p) (= tr^q^pq' 1 )) 

(where p £ T q (Sym£)is identified with a matrix £ Sym n ). 

- Show that Sym^ is isometric to the product x Y n , more precisely, one has 
an isometry: A £ Sym^ i— > (log det A, df ^ A ) £ (R, ncan) x Y„ (where can stands 
for the canonical metric ofR). 

- Show that q i — y q w arc isometry, that is, Y n as well as X n are Riemannian 
symmetric spaces. 

Exercise 2.3. For n = 1, Y n is R+ endowed with 

- Y2 is "a hyperbolic plane", i.e. a homogeneous simply connected surface with 
(constant) negative curvature. Compute this constant. Does this correspond to a 
classical model of the hyperbolic plane? 

- Show that the SL(n, M.)-action on Y n factors through a faithful action ofPSL(n, R). 

2.0. 2. Symmetric matrices vs quadratic forms. We guess it is worthwhile to seize the 
opportunity and clarify the relationship between quadratic forms and their represen- 
tations as matrices. Let E be a vector space, and as above Sym(E) and Sym + (E) 
its spaces of quadratic forms, and those which are positive definite, respectively. A 
basis (ej) of E yields a matricial representation isomorphism 

P £ Sym(E) ^p= (Pfa), P(e 3 )) i3 £ Sym n {= Sym(R n )) 

Justified by a latter use (see $8), the restriction to Sym + (E), will be denoted by: 
Q £ Sym + (E) i->- q £ Sym+. 

In fact, these representations depend only on the scalar product on E for which the 
given basis is orthonormal. 

Now, given (Q, P) £ Sym + (E) x Sym(E), that is a pair of a scalar product to- 
gether with a quadratic form on E, one associates a Q-autoadjoint endomorphism / : 
E — » E, representing P by means of Q, that is P(x, y) = Q(x, f(y)) = Q(f(x),y), 
where P and Q are understood here as symmetric bilinear forms. 

Fact 2.4. Given (Q, P) and their associates (q,p), the endormorphism f has a ma- 
trix representation A = q~ 1 p. 

Conversely, given Q and a Q-autoadjoint endomorphism f, its corresponding qua- 
dratic form P has a matrix representation p = qA. 

2.1. Flats. Let B = (e 4 ) be a basis of E. The flat F B C Sym + (E) is the space of 
scalar products on E for which B is orthogonal. It is parametrized by n positive reals 
Xi, its elements have the form: Ea^e* ® e*. 



4 



ABDELGHANI ZEGHIB 



Exercise 2.5. Show that the metric induced on F$ is given by Ej — and that 

X i 

(h,...,t n ) e (M n , can) i-> E exp(i;)e* ® e* G F B 

is an isometric immersion. 
Prove that: 

- Fb is totally geodesic in Sym + (E).(Hint: make use of the isometries of Sym^, 
p i — y o-*poi, where Oi is the reflection fixing all the ej, j ^ i, and <7j(e$) = — ej. 

- GL(E) acts transitively on the set of flats of the form Fq (Hint: relate this to the 
simultaneous diagonalization of quadratic forms). 

- Any geodesic of Sym + (E) is contained in some F&. 

3. An individual left invariant metric 
Standard references for this section and the following one are ll6ll9l [T7l[T8ll . 

3.0.1. Equation of Killing fileds. Let (M, (,)) be a pseudo-Riemannian manifold. 
A Killing field X is a vector field that generates a (local) flow of isometries. Any 
vector field X has a covariant derivative which is an endomorphism of TM defined 
by: D X X : u G T X M i-> V u X(x) G T X M, where V is the Levi-Civita connection 
of the metric. 

Fact 3.1. X is a Killing field, iff, D X X is skew-symmetric with respect to {, ) x , for 
any x G M. 

Proof. Let us first recall that this is the case for the Euclidean space: in fact, this 
is equivalent to that the Lie algebra of the orthogonal group is the space of skew- 
symmetric matrices. 

In the general case, assume x generic, that is, X(x) ^ 0, and consider N a small 
transverse submanifold (to X at x). Let Y and Z be two vector fields defined on N, 
and extend them on open neighborhood of x, by applying the flow of X, that is by 
definition: [X, Y] = [X, Z] = 0. 

If X is Killing, then (Y, Z) is constant along X: X.(Y, Z) = 0. Thus, by definition 
of the Levi-Civita connection, (VxY, Z) + (Y, V xZ) = 0. Apply the commutations 
V X Y = VyX, V X Z = V Z X, to get: = (VyX, Z) + (Y, V Z X), that is D X X is 
skew-symmetric. 

It is also easy to use those arguments backwards, that is, if D X X skew-symmetric 
for any x, then X is a Killing field. □ 

3.0.2. Three Killing fields. Let now X, Y and Z be three Killing fields. Apply 
skew-symmetry for all their covariant derivatives, and get (at the end of substitutions) 
the following formula: 



(2) 



2(V X Y,Z) = ([X,Y],Z) + (\Y,Z],X) - ([Z,X],Y) 
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Remark 3.2. Observe the beauty ( = symmetry) and easiness of the formula! 
Remark 3.3. Recall Koszul's formula for the Levi-Civita connection: 

2(V Y Z,X)= Y(Z,X) + Z(X,Y) - X(Y,Z)- 
(3) (Y,[Z,X]) + (Z,[X,Y]) + (X,[Y,Z}) 

It implies the previous formula of three Killing fields. Conversely, this last for- 
mula yields Koszul's one for any combination of Killing fields with coefficients (non- 
necessarily constant) functions on M. In particular (T2]) yields (O if M is homoge- 
neous. In fact, as (|2]) holds as X, Y and Z are pointwise Killing at order 1, © yields 
© also in the general non-homogeneous case. 

3.1. Left invariant metrics. We are interested now on Riemannian metrics on a Lie 
group G which are left invariant, i.e any left translation x \- > gx is isometric. This 
is in particular the case of any flow 4> l {x) = g f x, where {g t } is a one-parameter 

subgroup of G. Its infinitesimal generator X(x) = ^r(x)\t=o is a Killing field. This 
is a right invariant vector field: X(g) = X(l)g. 

Therefore, a left invariant metric is exactly a metric admitting the right invariant 
fields as Killing fields. 

Another characterization is that a left invariant metric is one for which left invari- 
ant fields have a constant length, (but they are not necessarily Killing). 

Such a metric is equivalent to giving a scalar product on the tangent space of one 
point in G, say T\G, i.e. the Lie algebra of G. We keep the same notation (, ) for 
both the metric on G and the scalar product on its Lie algebra Q. 

Here, to be precise, we define the Lie algebra Q as the space of right invariant 
vector field on G. (The other choice, i.e. that of left invariant vector fields, would 
induce modification of signs in some formulae). 

From the formula of three Killing fields, one sees that the connection on G is 
expressed by means of the scalar product on Q and the Lie bracket. In other words, 
one can forget the group G and see all things on the Lie algebra. For instance the 
Riemann curvature tensor is a 4-tensor on Q. The Ricci curvature is just a symmetric 
endomorphism of Q. 

3.2. The connection. Any Lie group has a canonical torsion free connection defined 
for right invariant vector fields by: 

V X Y = ~[X,Y}. 

Any other left invariant connection can be written VxY = \[X, Y] + C(X, Y), 
where C : Q x Q — )> C? is a bilinear map, and the connection is torsion free, iff, 
C(X, Y) = C(Y,X). In the case of the Levi-Civita connection of a metric, one 
deduces from the formula of three Killing fields 
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2(C(X,Y),Z) = ([X,Z],Y) + (X,[Y,Z]) 
(4) = (ad* x Y, Z) + {ad* Y X, Z) 

(Here as usually ad u v = [u, v], and ad* is its adjoint with respect to (, )). 
It then follows: 



(5) V X Y = ^([X,Y] + ad* x Y + ad* Y X 

3.2.1. Sectional curvature. The following formulae for various curvatures follow 
from Eq. (0]), see for instance (6l for detailed proofs. 

(R(X,Y)Y,X) = - l([X,Y},[X,Y])-^([X,[X,Y],Y) 

- l([Y,[Y,X},X]) + (C(X,Y),C(X,Y)) 

- {C(X,X),C(Y,Y)) + (Y,[[X,Y],X]) 

3.2.2. Ricci curvature. The tensor C disappears in the expression of the Ricci and 
scalar curvatures! 



Ric(X,X)= - -B(X,X)-([Z,X],X) 

(6) is 4 |[X,e 4 ]| 2 + is^([ ei , ej ],X) 2 
where, 

• B is the Killing form of Q, that is the bilinear form (X, Y) i-» tr{adxady), 

• (ej) is any orthonormal basis of (Q, (, )), 

• and finally, Z is the vector of Q defined by (Z,Y) = trace(ady), that is, it 
measures the unimodularity defect of Q by means of (, ). 

3.2.3. Scalar curvature. 

1 1 

(7) r = --Ejjlfe^ej]! 2 - -SjB(ei,ej) - \Z\ 2 

3.3. Warning: left vs right. A metric is bi-invariant i.e. invariant under both left 
and right translations of G if and only if its associated scalar product is invariant 
under the Ad representation. In this case, the connection is the canonical one given 
by V xY = \[X, Y], for X and Y right-invariant vector fields. As said above, this 
torsion free connection exists on any Lie group, but does not in general derive from a 
Riemannian or a pseudo-Riemannian metric (i.e. it is not a Levi Civita connection). 
Its geodesies trough the neutral element are one-parameter groups, and its curvature 
is given by R(X,Y)Z = \[[X,Y},Z], and has a Ricci curvature Ric(X,Y) = 
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\B{X,Y), where B is the Killing form (recall that a connection, not necessarily 
pseudo-Riemannian, has a Ricci curvature, Ric(X, Y) = tr(Z i-> R(X, Z)Y). 

3.3.1. Other quantities. The following fact, left as an exercise, gives characterization 
of bi-invariant metrics: 

Fact 3.4. A left invariant metric is right invariant (and hence bi-invariant) iff it sat- 
isfies one of the following conditions: 

(1) For any right invariant (Killing) field X, (X, X) is constant on G. 

(2) The orbits of any such X are geodesic 

(3) The orbit ofl(£ G) under any such X is a one parameter group. 

All this suggests the possibility to define other quantities essentially equivalent to 
the Ricci curvature. Say, in the bi-invariant case, the Ricci curvature is essentially the 
Killing form, and hence, in the general left invariant case, the remaining part Re (of 
the Ricci curvature) is an obstruction to the constancy of (X, X), or (equivalently) 
an obstruction for one parameter groups to be geodesic. A naive construction goes 
as follows. For X a right invariant vector field, let l x : x € G i— > (X(x), X(xj), its 
length function, and dlf its differential at I. Define Re(X, X) = tr{dl* dl* )... 



4. Curvature mappings on X n 

4.1. All scalar products together: The space Sym + (Q) et al. We are now consid- 
ering all left invariant Riemannian metrics on G. As said above the space of such 
metrics can be identified with Sym + (Q), the space of positive definite scalar prod- 
ucts on Q. Let Sym(Q) be the space of all quadratic forms on Q. Then, the above 
formula for the Ricci curvature determines a map: 

Ric : Sym + (Q) — > Sym(Q) 

4.2. Aut(Q) -action on Sym + (Q). Not only interior automorphisms of Q, but also 
"exterior" ones, i.e. general automorphisms act on Sym(Q). Their group Aut(Q) is 
sometimes identified with Aut(G), assuming implicitly that G is simply connected. 
Of course, this action is compatible with all mappings that will be discussed below. 

4.3. All Lie algebras together. We can now go a step further and deal with all Lie 
algebras of a given dimension n. As vector spaces they are identified with R n . The 
space of quadratic forms and the positive definite one are Sym n and Sym+. In short, 
for any Lie algebra Q (endowed with a basis allowing one to identify it with R n ), we 
have a Ricci and a scalar curvature mapping, involving the bracket structure of Q : 

Ricg : Sym^ — > Sym n 
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rg : Sym^ — > E 

4.4. Formulae. The Lie algebra Q has a basis (a). An element of Sym n is denoted 
by p = (xij)i<,i,j,<n- The structure constants Cfp are defined by [a, ef\ = C^eu- 

(Here and everywhere in this paper, if a letter is repeated as a lower and an upper 
index, like k is the last equation, we use the Einstein summation convention, for 
example, a\bj stands for Yija\bj, etc.). 

From Formula ©, we have: 

Ricg : q = (x ah ) G Sym+ H> p = (X ab ) E Sym n 

Where: 

(8) X ab = - X -B ah - l -C k a fi\ jXkl x^ + ^Cf k C] lXpa x qb x^x kl + X -E ah 

and, B a b = C 3 ai C l h - is the matrix representing the Killing form, (x^) is the inverse 
matrix of (xy), and 

E ab = Cl^ s (C l sa x lb + C l sb x la ) 

This last term depends on q, but vanishes identically if Q is unimodular. So, assuming 
Q unimodular will simplify and shorten the formula. 

4.5. Parameter. To simplify, we will restrict ourselves to unimodular algebras, and 
so (E ab ) disappears. Any Lie algebra is denned by a system (C^ ) which furthermore 
satisfies the Jacobi identity. We can then consider a mapping, 

Ric : 0, q) = ((C§), (x ab )) G M™ 3 x Sym+ ^p = (X ab ) E Sym n 

(9) X ab = Afifli - \c k ai C l b3 XM^ + \cl k C) [ x va x ab x^x kl 

This map is equivariant with respect to the GL(n, E) x GL(n, E) -action on the 
source and the GL(n, M)-action on the target. 

5. Three-dimensional case: Bianchi geometries (of class A) 

A Bianchi space (or geometry) is a homogeneous Riemannian 3-manifold to- 
gether with a transitive free action of a Lie group. This is therefore equivalent to 
giving a left invariant Riemannian metric on a 3-dimensional Lie group. These groups 
have been classified by Bianchi (see for instance [9, 20]). They split into classes A 
and B, according to they are unimodular or not. To simplify we will consider here 
only class A. 
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5.0. 1. Milnor (or Bracket cyclic) bases of a 3-dimensional Lie algebra. Let G( a ,b,c) 
be the Lie algebra generated by {u, v, w} with relations: 

[u, v] = aw, [v, w] = bu, [w, u] = cv 

It is easy to show that this is actually a Lie algebra, i.e. that the Jacobi identity is 
satisfied. 

Conversely, a basis B = {u, v, w} of a Lie algebra Q is called a Milnor basis 
(|9l) if it satisfies the previous relations: [u, v] = aw, [v, w] = bu, [w, u] = cv. In 
particular, Q is then isomorphic to Qub,c)- 

5.1. Invariance of Milnor flats under Ric. Let us recall here that in dimension 3, 
giving the Ricci curvature is equivalent to giving the full Riemann curvature tensor 
(in higher dimension, Ricci is too weaker than Riemann). 

Recall that the flat F& determined by B is the space of scalar products on Q for 
which B is orthogonal. 

We will say that F& is a Milnor flat if B is a Milnor basis. 

The flat Fg is parametrized by 3 positive reals x, y and z, where: 

x = (u, u) , y = (v, v) , and z = (w,w). 

The corresponding metric will be denoted by (x, y, z) G (R + ) 3 
Define the extended flat Fg to be the set of all quadratic forms which are diagonal 
in the basis B, i.e. they have the same form as elements of F&, but x, y and z are 
allowed to be any real numbers. 

One uses the orthonormal basis {^j, ^7=} for the metric (x, y, z), and com- 
putes from Formula ((6]) (or Q), 







Ric(u, 


u) = 


1\° yz 


"« 2 f 


-c 2 f) + ac 


(10) 


< 


Ric(v, 


v) = 


1 ( r 2 y 2 

2V C xz 




-6 2 f) + a6 






Ric(w, 


w) = 


2\ u xy 


-6 2 § 


-c 2 f) + 6c 



In fact, Ric is diagonal in the basis B, i.e. Ric(u, v ) = ... = 0. We can on the 
other hand perform some simplification, for instance, for Ric(u, u) we have: 

-(b 2 — -a 2 - -c 2 -) + ac = —(b 2 x 2 -c 2 y 2 -a 2 z 2 + 2acyz) 
2 yz y z 2yz 

= —(b 2 x 2 -(cy-az) 2 ) 
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Proposition 5.1. Consider a 3-Lie algebra Q, and a Milnorflat Fq C Sym^. Then, 
the Ricci map sends: 

(x,y,z) £F S 4 (X, Y, Z) G F~jg 

and is given by: 



X = 


2yz 


[b 2 x 2 — (cy 


-az) 2 ) 


Y = 


J_ 

2xz 


[c 2 y 2 - (az 


- bx) 2 ) 


Z = 


J_ 

2xy 


[a 2 z 2 - (bx 


- cy) 2 ) 



Remark 5.2. Observe the complete symmetry of these equations: there is for in- 
stance "a duality" x i->- b: everywhere the coefficient of x (resp. x 2 ) is b (resp. b 2 ). 
Recall that the coordinate x correponds to the vector u, and they are both related to the 
coefficient b, by the fact that the unique bracket proportional to u is [v, w] = bu. The 
same observation applies to the other variables, following a same correspondence: 
(x,y,z) h-> (b,c,a). 

5.2. Scalar curvature. Similarly, from Formula ©, we infer: 

(12) r = (—b x — c y — a z + 2acyz + 2abxz + 2bcxy). 

2xyz 

6. Structure of unimodular Lie algebras in dimension 3 

Proposition 6.1. Any unimodular 3-Lie algebra has a Milnor basis. More precisely 
(a, b, c) 6 M 3 i— > G( a ,b,c) ^ £ gives a parametrization of the space of unimodular 3- 
Lie algebras C. The diagonal action o/R* 3 on M 3 , and that of the permutation group 
S3 by permutation of coordinates, preserve the isomorphism classes of algebras. 

In fact, this is equivalent to the more precise: 

Corollary 6.2. There is six isomorphism classes of unimodular 3-algebras, repre- 
sented as follows, where below + (resp. —) means any positive (resp. negative) 
number, e.g. +1 (resp. —1). 

(1) (0, 0, 0): the abelian algebra R 3 . 

(2) (0, 0, +); the Lie algebra of G = Heis, the Heisenberg group. 

(3) (0,—,+): the Lie algebra of G = Euc, the group of rigid motions of the 
Euclidean plane (i.e. the isometry group of the affine Euclidean plane). It is 
a semi-direct product R x R 2 , where R acts on R 2 , by 

I cos t — sin t \ 
V smt cost J 
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(4) (0, +, +)• the algebra of the group G = SOL, the group of rigid motions 
of the Minkowski plane (i.e. the Minkowski space of dimension 1 + 1), G = 
R ix R 2 , where R acts on R 2 by 



They are labeled respectively, Bianchi: I, II, VIIq, VIq, VIII and IX. 

Proof, (see ifTTl ). In dimension 3, there are exactly two semi-simple Lie algebras, 
so(3) and sl(2, R), and there are no semi-simple Lie algebra of dimension < 2. 
Therefore, if a 3-algebra Q is not semisimple, then it contains no semisimple algebra, 
that is, Q is solvable. 

Assume Q contains an (abelian) ideal isomorphic to R 2 . Take any supplementary 
one dimensional subspace (hence a subalgebra) R. Then, Q writes as a semi-direct 
product of R acting on R 2 . Since Q is assumed to be unimodular, this action is via 
a one parameter group of SX(2,R). These one parameter groups split into: elliptic, 
parabolic and hyperbolic types. We obtain respectively: Euc, Heis and SOL. 

It suffices therefore to show the existence of such an ideal R 2 . For this, let us 
remark that there exists always an abelian ideal of dimension 1, say A = R. Indeed, 
if the commutator ideal has dimension 2, then since it is solvable, its commutator has 
dimension 1 or 0... 

Now, since Q acts on A, the Kernel C has at least dimension 2 and contains A. 
Let us consider here the case dim£ = 2, since the dimension 3 case is easier. By 
definition (of the Kernel) this is an ideal, and since it has dimension 2 and has a 
non-trivial center (it contains A), then it is abelian. □ 



In this section, we present a general setup where the previous constructions can be 
defined. In other words, we "summarize" how one associates to a Lie group various 
rational dynamical systems. 

7.1. A rational map. ( £14.41) . Let Q be a Lie algebra of dimension n, with a basis 
(ej), such that [e$, ef\ = C^e fc , and assume to simplify that it is unimodular. Then, 
we have a rational map: 




(5) (-,+,+): aZ(2,R). 

(6) (+,+,+): so(3). 



7. Summary; further comments 



Ricq : (x ah ) e Sym+ H> (X ab ) <E Sym, 



(13) X a b — — — Bab — —CaiCljXkiX 1 -* + —CffcC^XpaXqbX^X' 



kl 



as well as a rational (scalar curvature) function: 



(14) 



rg : {x ab ) G Sym+ H> X ab x 



12 



ABDELGHANI ZEGHIB 



(recall that (x 1 - 7 ) is the inverse matrix of (xij)). 

7.2. Aut(Q)-equi\ariance. Both Ricg and rg are respectively equivariant and in- 
variant under the Aut{Q) -action on Sym n (identified with Sym{Q), the space of 
quadratic forms on Q). 

7.3. Extensions. Actually, everything extends to left invariant pseudo-Riemannian 
metrics on G, or equivalently to Sym*(Q) the space of non-degenerate quadratic 
forms on Q. The same formulae allow one to calculate the Ricci and scalar curvature 
of such metrics. 

• Now, the formulae may have sense even for some degenerate quadratic forms ! 

• We can also consider complex quadratic forms on the complexification of Q. 
They form a space Sym(Q) <8> C identified with Sym n {£), the space of symmetric 
complex matrices. In other words, to a Lie group G of dimension n is associated 
an Aut(G) -equivariant rational transformation Rico on Sym n (C). We also have an 
^^(GQ-invariant meromorphic function tq- 

• Ric G can be written as a rational vectorial map (^)j : c n(n+1)/2 -)• c n(n+1)/2 
where Pj and Qi are homogeneous polynomials of a same degree 2n (on n(n + l)/2 
variables). 

• From this, one gets a rational transformation of the projective space cp n ( n + 1 )/ 2 -' 1 
PSym n (C). We will denote it by Ric^, to emphasize that it is the extension to the 
complex projective space. For instance, if n = 3, then we have a rational map on 

CP 5 . 

• Rico is equivariant under scalar multiplication. In fact, as this follows from 
its defining formula, Rico is polynomial when restricted to matrices with det = 1. 
Therefore, a representative of Ric^ (on the projective space) is the polynomial map 
(Pj)j (of degree 2n). 

• Rice is invariant under the (algebraic) action of (the complexification of) Aut(Q) 
on cP n ( n+1 )/ 2_1 (= PSym(Q)®<C). It then determines a map on the quotient space. 
It depends however on the meaning to give to such a quotient space (by Aut(Q)). As 
an algebraic action, it has a poor dynamics, and a nice quotient space can be thus con- 
structed for it. There is in particular a notion of "algebraic quotient". In dimension 
n = 3, the algebraic quotient space has dimension 5 — 3, more exactly, it is a (singu- 
lar) compact complex surface Sq, say. We have then associated to a 3-dimensional 
Lie group a rational map on a compact complex surface Sq- This map seems to have 
a "poor dynamics", for instance, it has in general a vanishing entropy. We guess 
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nevertheless that (other) "dynamical invariants" of it can characterize the group G 
(i.e. two different groups have different invariants). It is also worthwhile to see what 
happens in higher dimension case. 

7.4. Forget invariance. In the formula defining Rice we can consider any system 
of parameters (C^), not necessarily satisfying the Jacobi identity of Lie algebras. We 
obtain a big family of rational transformations generalizing those associated to Lie 
groups. In this case, various dynamical types may appear. We think it is worthwhile 
to investigate the structure of this parameter space, and to understand inside it, the 
(algebraic) set of Lie algebras, the algebraic actions on it... 

7.5. Cross sections, Flats. Let us call a cross section S for Ricq or Rice a sub- 
manifold in Sym+ (resp. PSym n (C)) which is invariant under Rice (resp. Ric^) 
and such that S meets any G-orbit in a non-empty discrete set. The last condition im- 
plies in particular that S is transversal (at least in a topological sense) to the G-orbits. 
In fact there are weaker variants of this definition which can be useful, in particular 
in a geometric algebraic context. 

• It is in the case where G is unimodular and has dimension 3 that cross sections 
occur easily. In this case, one can in fact find them, as affine (resp. projective) 
subspaces for Rice (resp. Ric^) of dimension 3 (resp. 2). Indeed, as explained in 
$61 the Lie algebra of such a group has a Milnor basis B = {u, v, w} ( §5.0. 1 b . and the 
flat Fb (£ i2.1l) . or more formally its "extension" Fg, i.e. the space of quadratic forms 
diagonalizable in £>, is invariant under Ricq- In order to see that one gets in this way 
a cross section, it remains to show the abundance of Milnor bases as in the following 
statement, 

Exercise 7.1. Prove that any quadratic form can be digonalized in some Milnor basis 
of Q. (Hint: this can be done by checking case by case. For instance, for the group 
SOL, its Lie algebra is generated by X, Y, Z, with relations [X, Y] = Y, [X, Z] = 
—Z and [Y, Z] = 0. Consider u = X+T, where T belongs the the plane V generated 
by Y and Z. Then, the restriction of ad u on V, satisfies ad\ = —1. Choose u to be 
orthogonal to V (with respect to the given metric). Consider a non-vanishing vector 
v G V, and let w = ad u (v). Then {u, v,w} is a Milnor basis, because of the fact 
ad\ = 1. We claim that v can be chosen such that w is orthogonal to v. This is a 
calculation in the basis {Y, Z}). 

• Let us point out the following polynomial presentation of Ricq. As was said 
above, Ricq is invariant under scalar multiplication, that is, it suffices to consider 
its restriction on unimodular matrices SSym n , in which case, it becomes polyno- 
mial. In particular, from £ j5.1[ Ric^ has the following form as a cubic homogeneous 
polynomial map: 

(15) Ric G : C 3 -> C 3 
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(a?, y, z) i y -(x(b 2 x 2 — (cy — az) 2 ), y(c 2 y 2 — (az — bx) 2 ), z(a 2 z 2 — (bx — cy) 2 )). 

7.6. Formula in each case. 

7.6.1. Case of SO(3): a = b = c=l. 
Ric(x, y, z) = ^(x(x 2 - (y - z) 2 ),y(y 2 - (z - x) 2 ), z(z 2 

7.6.2. Case of SL(2,R): a = b = 1, and c = -1. 

Ric(x, y, z) = ^(x(x 2 - (y + z) 2 ),y(y 2 - (z - x) 2 ), z(z 2 

7.6.3. Case of the Heisenberg group Heis: a = b = 0,c = 1 

I 

Ric(x,y,z) = -y 2 (-x,y,-z). 

7.6.4. Case of Euc: a = 0,b = —1, c = 1 

Ric(x, y, z) = ^(x(x 2 - y 2 ),y(y 2 - x 2 ), -z(x + y) 2 ) 

7.6.5. Case of SOL: a = 0, b = c = 1 

Ric(x, y, z) = ^(x(x 2 - y 2 ),y(y 2 - x 2 ), -z(x - y) 2 ) 

1.1. Bianchi-Ricci flow. Recall that the Ricci flow associated to a compact manifold 
M (of finite volume) is an evolution equation on its space Met(M) of Riemannian 
metrics: 

at n 
where n = dimM and < r(g) >= f r(g)dv g /Vol(M,g) is the average scalar 
curvature of q ||9l. 

7.7.1. The vector field IZicc- Now, if G is an n-dimensional Lie group, then this 
gives a classical differential equation on the space of its left invariant Riemannian 
metrics, where one takes a punctual value of the scalar curvature instead of its average 
(since this scalar curvature is constant). Equivalently, this is a vector field on SymX. 
In fact, all this is derived from our previous rational map Rice- Since Sym^ is an 
open set in the vector space Sym n , the vectorial map Ricq : Sym^ — > Sym n can 
be alternatively seen as a vector field on 5ym+, say, TUcq. 

Observe that IZicc is invariant under the ylut(G)-action on Sym+ (which is 
equivalent to the fact that Ricq is equivariant under the (linear) action of Aut(G)). 



-(x-y?)). 

-(x + y) 2 )). 
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Let us denote a generic point of Sym^ by q, and consider the radial vector field 
V(q) = q. The previous differential equation, which we will call the Bianchi-Ricci 
flow associated to G is the vector field 

-2Kic G + 2— V 
n 

7.7.2. Commutation. Consider the bracket [IZicc, V] = DyRicc — Djh cg V, where 
D u denotes the usual derivation in the u-direction. This equals — Rice, since, 
DyRicc = 0, i.e. the map Rice is invariant under multiplication; and DV = 
Identity, everywhere. Therefore, [TZicG,V] = —TUcq. Because of this commu- 
tation rule (that is, the two vector fields generate a local action of the affme group), 
the essential dynamics of the Bianchi-Ricci flow comes from the TZicc-psit- 

7.7.3. Bianchi-Hilbert-Ricci flow. The remark applies to any combination of IZicc 
and V: understanding one combination allows one to understand the others. A famous 
one is Sin = TZicc — §V, which can be called in this context the "Bianchi-Einstein 
flow", since the tensor Ric(g) — |g of a Riemannian manifold (M, g) is called Ein- 
stein tensor (this is, essentially, the unique combination of Ric(g) and g which is 
divergence free). However, in order to prevent confusion with "Einstein equations" 
and some related flows which will be considered below, Sin could be better called 
Bianchi-Hibert flow. Indeed, the function 

% : p G Sym^ i-> r(p)y / det(p) G M 

is the substitute of the classical Hilbert action in the case of left invariant metrics. 
Indeed: 

Exercise 7.2. Show that Sin is a gradient vector field. More exactly, Sin = VH, 
where the gradient V is taken with respect to the metric of Sym^. 

Remark 7.3. The computation can be handled in a more explicit way on a Milnor 
flat Fq ( £15.11 ). where the Hilbert action has the form: 

H(x, y, z) = — ■ (—b 2 x 2 — c 2 y 2 — a 2 z 2 + 2acyz + 2abxz + 2bcxy) 

2^/xyz 

and the metric is 

dx 2 dy 2 dz 2 
x 2 y 2 z 2 

7.7.4. Restriction on SSym^- The interest of the normalization in the definition of 
the Ricci flow is to let it preserving the volume of the Riemannian metric, that is, the 
total volume remains constant under evolution. In the case of left invariant metrics, 
this is equivalent to the fact that the vector field TZicc — is tangent to SSym n . 

This in turn is equivalent to the fact, that for any q £ SSym n , Ricc{q) — 

is trace free, which follows from the very definition of the scalar curvature r. This 

allows one to justify the following simplification: write equations assuming q G 
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SSym n , i.e. det(g) = 1, which gives polynomial equation. However, in order 
to keep this polynomial natural, do not take reduction of variables from the equation 
det(g) = 1. To be more concrete, consider a flat Fg, then instead of the rational forms 
of Ricc{x,y, z) and rc(x,y,z), we assume xyz = 1 which leads to polynomial 
forms: § §5. II and Formula [T2l (but we do not go further and eliminate one variable, 
say z = -M. We can then write the Bianchi-Ricci flow as follows 

xy 

(16) 

x(^bx(-2bx + cy + az) + \(?y 2 + ^a 2 z 2 - \acyz) 
y(^cy(-2cy + bx + az) + \b 2 x 2 + \a 2 z 2 - \bcxz) 
z(^az(-2az + bx + cy) + \b 2 x 2 + \(?y 2 - \bcxy) 

7.7.5. Differential equations on a projective space. In the same way, we associate 
to a Lie group G an Aut(G)-invariant one dimensional complex algebraic foliation 
on the projective space PSym n (C). Here, among combinations of the vector fields 
IZicc and V, only IZicc is relevant, since the radial vector field V becomes trivial 
on the projective space. In the case of a unimodular 3-group, we have the following 
homogeneous cubic differential system on C 3 : 

' f t = x (b 2 x 2 - (cy - az) 2 ) 

( 17) S= v((cV-(«-fa0 2 ) 

£| = z(aV - (ta - cy) 2 ) 

i. 

7.7.6. Dynamics, compactifications. It is the dynamics of the Bianchi-Ricci flow 
—21Zicc + ^-V which was investigated in the literature (9j[13l. As we argued above, 
this is essentially the same as that of the Einstein-Hilbert field VH. But, as a gradient 
flow, its dynamics is completely trivial on Sym^... The point is to study the behavior 
of orbits when they go to an infinity boundary doo Sym+. There is however several 
ways to attach such a boundary to (the non-positively curved Riemannian symmetric 
space) Sym£. One naturally wants to interpret ideal points as collapsed Riemann- 
ian metrics. With respect to this, the Hadamard compactification seems to be the 
most pertinent (see for instance ifToll ). On the other hand, the advantage of algebraic 
compactifications (e.g. the projective space) is to extend the dynamics... 



-21Zic G +2 r -V{x,y } z) 
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8. Hamiltonian dynamics on Sym+ 

After consideration of some maps and vector fields, we are now going to study 
second order differential equations on Sym£, the prototype of which is the geodesic 
flow of Sym+, and then the "Einstein flow" associated to a Lie group. 

8.1. Geodesic flow. Write the metric on Sym+ as: L(q,p) = (p,p) q = tr(q~ 1 pq~ 1 p). 
Since Sym^ is open in Sym n , its tangent bundle trivializes TSym+ = Sym+ x 
Sym n . We will use the usual notations |^ for the horizontal and vertical differ- 
entials d p L and d q L. 
We have: 

dL 

_(<5g) = tr(-q' l (5q)q~ 1 pq' l p-q' 1 pq' 1 {5q)q' l p) = -2tr{(5q)q' 1 pq' 1 pq' 1 ) 
where 5q is a horizontal tangent vector, i.e. an element of Sym n . 

r)T 

-^(5p)=2tr((5p)q~ 1 pq- 1 ). 

Now, write: q = q(t), p(t) = q = and compute 

^~{&P) = 2tr((5p)q- 1 [-2 q q- 1 q + q]q- 1 ). 

The Euler-Lagrange equation is obtained by taking Sq = 5p = A, and writing for 
any A, 

dt dp dq 

This reads: 

tr{2Aq- l {-qq- l q + q)q~ l ) = 0, \/A G Sym n . 

and therefore, 

Fact 8.1. The equation of geodesies of Sym+ is the second order matricial equation 
on Symf L : 



q = qq 1 q 



or equivalently (in the phase space): 



(18) 



P 



p 

pq~ 1 p. 
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8.2. Other pseudo-Riemannian and Finsler metrics on Sym+. There is a canon- 
ical GL(n, R)-invariant form uj on Sym+: 

co q (p) = tr{q~ l p). 

We can then associate to any reals a and (5 a Lagrangian: 

L a ,/3(q,p) = a(uj q (p)) 2 + P(p,p) q = a{tr{q~ l p)f + ptr^pq^p). 

For generic a and (3, this is a homogeneous pseudo-Riemannian metric, but it can 
degenerate for some values. 

Similarly, there are homogeneous Finsler metrics: 

Fa,p{q,p) = au q {p) + fiyj (p,p) q = a^r^p)) + fi^tr(q- l pq- l p. 

Exercise 8.2. Write the Euler-Lagrange equation for L a>( g and F a ^. 
- Solve the geodesic equation for Symi- 



9. Einstein Equations in a Gauss gauge 



Cylinders. Let M be a differentiable n-manifold endowed with a family of Rie- 
mannian metrics g t , t is a "time" parameter lying in an interval I. Consider the 
Lorentz manifold M = I x M endowed with the metric 

(,) =9 = ~dt 2 + g t , i.e. g~( t ,x) = ~dt 2 + (g t ) x 

Such a structure is sometimes called a cylinder. Our purpose is to relate geometric 
(e.g. curvature) quantities on M and M. For a fixed point (t, x), R, Ric, and r will 
denote the Riemann, Ricci and scalar curvatures of (M, g t ) at x and V its Levi-Civita 
connection. The corresponding quantities for M are noted by R, Ric and f and V. 

9.1. Second fundamental form. The (scalar) second fundamental form of {t} x M 
is denoted k t (or sometimes simply k). Actually, the second fundamental form is 
defined as a vectorial form: II(X, Y) equals the orthogonal projection of V xY on 
Reo, where eo = Jj. 

The scalar second fundamental form is defined by 

k(X,Y) = (II(X,Y),e ) = (V x Y,e ). 

The Weingarten map a = a eo is defined by: 

a(X) = -V x e . 
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We have: 

k(X, Y) = (V X Y, eo) = X(Y, e ) - (V x e , Y) = + (a(X), Y). 

In other words, a is the symmetric endomorphism associated to k by means of 
the metric g (we will use sometimes the notation at as well as g t and k t , in order to 
emphasize the dependence on i).(The definition of k and a coincides with that in the 
Riemannian case. The unique difference is that here, II = —/ceo, since eo is unit 
timelike, i.e. (eo,eo) = — 1). 

9.2. Geometry of the product. Consider e 1; . . . , e n a frame of vector fields on M, 
that we also consider as horizontal vector fields on M. By definition, they commute 
with e (= Jj). 

Fact 9.1. We have: 

(19) V eo e 

(20) kt 

(21) (R(e ,ei)ei,e } 

(22) Ric(e ,e ) 



= (the trajectories ofeo are geodesic), 

= Qj_(at{ei),ei} + (at(ei),ei), 
= i-Ma t ) + tr(a 2 t). 



Proof. • We have 

= d/dt(e ,ei) = (V eo e ,ej) + (e ,V eo ei). 

But 

(e , V eo ei) = (l/2)e;.(eo,e ) = 0, 
since eo and e, commute. Therefore (V eo eo, ei) = 0, Vi 

• We have 

d 

Q^9t(ei,ej) = e Q {e j ,e j ) = (V eo ei,ej) + (V eo ej,ej) 
Since eo commutes with e^ and ej, this also equals: 

-((a(ei),ej) - (a(ej),ei) = -2k t {ei,ej) 

• Computation of R(eo, ej)ej, eo): Because of the commutation relations, and be- 
cause eo is geodesic, we have, by definition of the curvature: 

R(e ,ei)e = -V eo a(ej), 

and thus: 

(R(eo,ei)e ,ei} = -(V eo a(ej),ej) = -e (a(ej), e^) - (a(ej), a(ej)), 
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and since a is symmetric, this also equals: 

d 

-7^(a(e,),ej) - (a 2 (e;),ej) 

And hence, 

d 

(R(e ,ei)ei,eo) = — (a^e^e;) + (a 2 (ej),ej) 

• We can assume that at a fixed point (t, x), the basis (ej)i>i is orthonormal, and 
taking the sum (over i > 0) we get: 

d 

Ric(e , e ) = ^ r (°t) + tr ( a t )• 

□ 

Remark 9.2. In fact, the meaning of "Gauss gauge" is nothing but that e is unit and 
has geodesic orbits. 

9.3. Gauss equation. It describes the relationship between the sectional curvatures 
for R and R: 

(-R(cj) Cj)ej, &i) = (i?(ej, j Gi) + k{ei, ei)k{ej , ej) 
(23) -k(ei,ej)k(ei,ej) 

(observe this difference of sign of the /c-term, in comparison with the Riemannian 
case). 

9.4. Einstein evolution equation for k t . Again, assume (ej) orthonormal, fix i, and 
take the sum over j > 0. We first have: 

T,jk(ei, ei)k(ej,ej) = fc(ej, ei)tr(a) = ir(a)((a(ej), e^)) 

and 

S i /c(e i ,ej)/c(e i ,ej) = (a 2 (ej),ej)) 

(Indeed in matricial notations, = ajj = k(ei,ej), and thus (a 2 )„ = ^jCiijCiji). 
Therefore, if we consider the quadratic form /, defined by: 

^(ci)Cj) = Sj (/c(ej, ej)/c(ej , ej) — /c(ej, ej)/c(ej, ej)) 

then its associated endomorphism is: 

tr(a)a — a 2 

• Ric(ei,ei) equals the trace of u — > R(u,ei)ei. Remember, (e^) is a Lorentz 
orthonormal basis, i.e. (ej, ej) = 0, for i / j, (eo, eo) = — 1 and (ej, ej) = +1, for 
j > 0. It then follows that 

R{ei,ei) = ^( J R(ej,ej)ej,ej) - (^(eo, e,)ej, e ) 

j>0 
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• Returning to the Gauss equation (I23I ). and taking the sum over j > 0, we get: 

Ric(ei,ei) + (R(e , e^e*, e ) = Ric(e^e^) + ((tr(a)a - a 2 )(ej),ej) 

• Replacing {R(eo, ej)eo, e«) by its previous value: 

Ric(ei,ei) + — (a(ei),ej) + (a 2 (ej),ej) = Ric(ei,ei) + ((tr(a)a - a 2 )(ej),ej) 
Equivalently, for any I,7e TM: 

Y) = -ffic(X, Y) + Ric(X, Y) + ((ir(a t K - 2a 2 )(X), Y>. 

Fact 9.3. Define the square power kt gt kt to be the quadratic form associated by 
means of gt with the matrix a 2 (where at is the matrix associated to kt via gt). Then: 

d_ 

dt 

9.5. Gauss constraints. Consider again the Gauss equation and take the sum over 

i,j > 0: 

f - 2Ric(e ,e ) = r + (tr gt k t ) 2 - tr(k tgt k t ) = r + (tr gt k t ) 2 - \k t \ 2 gt . 

9.6. Matrix equations. We are now going to write equations by means of symmet- 
ric matrices associated to the quadratic forms g t and k t (see §2.0.2b . For this, we fix 
x and a time to and choose an orthonormal basis (ej(to)) of T X M. We denote by q t 
(resp. p t ) the matrix associated with g t (resp. —2k t ), and by rict and rict (or simply 
ric and rid) those associated with Ric and Ric (recall they are the Ricci curvatures 
of respectively M, at (t, x), and (M, g t ), at x). With this, we have: 

• Evolution equations: 



-k t = -Ric + Ric + tr gt (k t )k t - 2k tgt k t . 



(24) 



q= p 

p= —ric + ric+\tr(q~ 1 p)p—\q' 1 pq~ l p 



• Gauss constraints (actually said Hamiltonian constraints) 

r — 2ric(eo, eo) = r + (£?"(<? -1 p)) 2 — tr(q~ 1 pq~ 1 p) 
= r + tr(q~ l pf - {p,p) q 
= r-L_i,i(g,p) 

where L_i i is the pseudo-Riemannian metric defined in JJU for the value (a, /?) = 

(-1,1) 
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10. BlANCHI COSMOLOGY 

We will now restrict ourselves to the vacuum case, i.e. M is Ricci-flat: Ric = 0, 
and thus also f = 0. We will also assume M is a Lie group G and the metrics on it (i.e. 
gt) are left invariant. Therefore such a metric is identified with an element q G Sym^ 
(n = dimG, the identification of Sym + (Q) with Sym+ comes from a choice of a 
basis). Now, ric becomes a map ric : Sym^ h-> Sym n , and r : Sym n \-> E. We get 
the (beautiful) ODE system with constraints: 



(25) 



V 



, L-i,i(q,p) 



ric(q) + \tr(q 1 p)p — \q 1 pq 1 p 
r(q) (Hamiltonian constraint) 



Remark 10.1. Observe that ric and r are basic functions, they depend only on q 
(and not on p). 



10.1. Isometric G-action on M. Here M = I X G, with g = —dt 2 + g t . A left 
translation x G G \- > hx is isometric for all the metrics g t , and therefore is isometric 
for g as well. 



10.2. The Bianchi-Einstein flow along and on a flat. Actually, there are other con- 
straints to add to the ODE system above, in order to get what we will call the Bianchi- 
Einstein flow. These (momentum) constraints will be considered below. Before, let 
us consider a subsystem of it, the restriction (of everything) to a Milnor flat Fq. The 
following proposition derives from Formulae (TTTT t and (fT2l) . 



Proposition 10.2. Let Fg be a Milnor flat, and TFq its tangent bundle, a point of 
which is denoted by (q,p), q = (x, y, z), p = (x 1 ', y' , z'). The Bianchi-Einstein flow 
on Fq is the following system of ODE on TFq, together with one algebraic constraint 
defined by a Lorentz metric on TFq and a basic function on Fq. (The phase space 
has thus dimension 5, and is a fiber bundle over Fq): 
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(26) 



x 



2±j(6V - (cy 



2xz\ L y 



az 



4\ x 



lx" 
2~P 



az-bxf) + \(*- + y- + *i)y>-\ 



2~p 



^ y (a 2 z 2 -(bx- C y) 2 ) + ^- + ^ + ^z' 



lz" 
2l? 



The phase space is a hypersurface ( maybe singular) N ( in TFq ) defined by the 
Hamiltonian equation: 



(27) l(x, yi z)(x',y',z') 
Where I is the Lorentz. metric ( on F$ ): 

(28) 



r(x,y,z) 



7 I i i i\ x y x z 

l (x,y,z){X ,V ,Z) = 1 

V >»> I X y XZ 



+ 



yz ' 



and r is given by Formula (T2v : 
1 



r(x,y,z) 



2xyz 



b 2 x 2 — c 2 y 2 — a 2 z 2 + 2acyz + 2abxz + 2bcxy) 



Exercise 10.3. Show explicitly that the constraint is preserved by the dynamics, i.e. 
the vector field determined by the differential equations is tangent to the "submani- 
fold" N C TFq defined by the constraint. 

10.3. Codazzi (or Momentum) constraints. The Codazzi equation establishes a 
relation between the intrinsic and extrinsic curvatures of a submanifold M in a Rie- 
manniann manifold M, and is in fact valid in the general background of pseudo- 
Riemannian manifolds provided the induced metric on the submanifold is also pseudo- 
Riemannian, i.e. it is not degenerate. More precisely, it states that some "partial 
symmetrisation" of the covariant derivative of the second fundamental form (all this 
depends only upon data on M) equals the normal part of the Riemann curvature ten- 
sor (this depends on M). The equation gives obstructions for a (vectorial) 2-tensor to 
be the second fundamental form of a submanifold. 
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In the case where M is a CMC spacelike hypersurface (i.e. with a constant mean 
curvature) in a Ricci flat Lorentz manifold M, one can deduce from Codazzi equa- 
tion, by taking a trace, that the second fundamental form k is a divergence free 2- 
tensor. This applies in particular to our case: our hypersurfaces are G-orbits and thus 
are CMC. 

Let us recall some definitions. Firstly, if A; is a symmetric 2-tensor on M, then its 
covariant derivative Vxk with respect to a vector X, is a 2-tensor: 

(V x k)(Y, Z) = Xk(Y, Z) - k(X7 x Y, Z) - k(Y, X7 X Z) 

Now divk is a 1-form, the trace of Vk (with respect to the metric of M), i.e. if 
(ej) is an orthonormal basis: 

divk(X) = ZiV e Mei,X) 

10.3.1. Divergence of left invariant quadratic forms on Lie groups. At first glance 
one can guess that left invariant objects are divergence free (with respect to left in- 
variant Riemannian metrics). This is however false (apart from some trivial cases). 

Let G be a 3-dimensional unimodular Lie group, endowed with a left invariant 
metric (, } = q G Sym + (Q), with a Milnor g-orthonormal basis {u, v, w}: [u, v] = 
aw, [v,w] = bu and [w,u] = cv (see £ j5.0.1l ). The proof of the following facts and 
corollaries is left as exercise. Let p G Sym(Q) represent a left invariant quadratic 
form. 

Fact 10.4. Let X, Y and Z be right invariant vector fields, with X(l) = e £ Q. Let 
g l = exp te. Then the derivative X.p(Y, Z) at 1 £ G is given by: 

X.p(X, Y) = ^- t p(Adg t (Y),Adg t (Z)) = p([X, Y],Z)+ p(Y, [X, Z)) 

Fact 10.5. For the basis {u, v, w}, we have: 

• V u u = V v v = V w w = 0, 

• 2V u w = (— c + a — b)v, 2V v w = (b — a + c)u ... (Use Formula (0) 

Corollary 10.6. Consider the left invariant quadratic form, pyi = du®dv + dv®du. 

Then: 

• 

u.pi 2 (u,e) + v.pi 2 {v,e) + w.pi 2 (w,e) = -c-a,fore = w 

= 0, for e = u, or e = v 

• Pii{u, y u w) + pi2(v, V v w) + pn(w, V w w) = 

• It then follows that uj = divpi2 is such that uo(u) = oj(v) = 0, and u>(w) = 
— (c + a), that is uj = — (c + a)dw. 
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Corollary 10.7. Let us say a Milnor basis is generic if (a + c)(a + b)(b + c) 7^ 0. 
Then, for a generic Milnor basis, any divergence free left invariant quadratic form 
(with respect to the metric for which this basis is orthonormal) is diagonalizable in 
this basis. 

In other words (keeping the previous notation), along a Milnor flat Fq, an element 
p G TqSyrnt satisfies the momentum constraints, iff, p G T q (Fg) (or in more linear 
words, p G Fq) 

10.4. Cross sections for the Bianchi-Einstein flow. On TSym^, the group Aut(G) 
acts, preserving the Bianchi-Einstein flow (determined by G). A cross section ( ^7.51 ) 
will play the role of a flow on a quotient space (for the Aut(G)-action). 

Proposition 10.8. The Bianchi-Einstein flow on a generic Milnor flat is a cross sec- 
tion of the full Bianchi-Einstein flow (with constraints) on Sym^ endowed with the 
Aut(G)-action. Generic flats exist except in the abelian and nilpotent cases, i.e. 
when G is R 3 or the Heisenberg group Heis. 

Proof. Firstly, one easily sees that if there exists a Milnor basis for which a and b 
/ 0, then after re-scaling, this basis becomes generic. This exists exactly when G is 
different from R 3 , and Heis. If we are not in these cases, then we can assume, after 
re-scaling if necessary, that all Milnor bases are generic. Let B be such a basis, and 
(q>p) G TSym+, then up to application of an element of Aut(Q), q G Fg. But since 
B is generic, if p satisfies the momentum constraints, then p G TFq, which means 
that TFg is a cross section. □ 

10.4.1. Case of G = R 3 . In the case of the Heisenberg group, there is exactly one 
momentum constraint which gives rise to invariant sets of the system. There is no 
such constraint in the case of R 3 , where we obtain the following system: 

q= p 

(29) { p = \tr(q~ 1 p)p — \q~ X pq~ X p 

L-i,i(q,p) = (the lightlike cone bundle of 
The spacetime M has a metric 

g = -dt 2 + t 2pi du 2 + t 2p2 dv 2 + t 2pS dv 2 . 

This is called a Kasner spacetime (observe that in some cases, e.g. pi = P2 = P3, 
this is just the Minkowski space) 11211711251. 



Exercise 10.9. Prove the previous form of g and solve the same problem in the case 
of the Heisenberg group. 
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10.5. Isometry group of M. As said previously, the left action of G on itself in- 
duces, by definition of its metric, an isometric action on M. In fact, if for some level 
(M, gt), there are extra-isometries (i.e. other than left translations), then they extend 
to M. More precisely, if the metric at some level, say t = 0, is identified with qg £ 
Sym + (Q), and M corresponds to a point (qo,Po) £ TSym^ = Sym+ x Sym n , 
and K C Aut{Q) is the stabilizer of (qo,Po), then, on the one hand, K acts as an 
isometric isotropy group for (M, go) (go corresponds to qo). On the other hand, K 
preserves the Bianchi-Einstein trajectory of (qo,po), and thus acts isometrically on 
M (as isotropy for any point identified with 1 £ G). 



10.6. An example: Bianchi IX. This means G = 50(3), or more precisely its 
universal cover the sphere 5 3 . In this case, there are Milnor bases with a = b = c = 
1. Any other Milnor basis satisfies these equalities, up to re-scaling. Also, all such 
bases are equivalent up to conjugacy and re-scaling. Yet, this is the most challenging 
case of Bianchi cosmologies (see for instance EH). As an example, TAUB-NUT 
spacetimes are exact solutions of the Bianchi-Einstein equations of class IX. They 
are characterized among Bianchi IX spacetimes as those having extra-symmetries, 
i.e. a non-trivial isotropy, which then must be 50(2) (and thus their isometry group is 
5 3 x 50(2), up to a finite index). Nevertheless, their high complexity (at least among 
exact solutions) led people to describe them as "counter-examples to everything"! In 
a Milnor flat where a = b = c = 1, these spacetimes correspond (up to isometry) to 
x = y, and x' = y'. The left invariant metric on G (at any time) corresponds to a 
Berger sphere, i.e. (up to isometry) a metric on the sphere derived from the canonical 
one, by rescaling the length along the fibers of a Hopf fibration. In other words, the 
set of solutions of the Bianchi-Einstein flow, which are Berger spheres at any time, is 
closed and invariant, say, the TAUB-NUT set. 



10.7. Effect of a non-vanishing cosmological constant. Instead of requiring M to 
be Ricci-flat, let us merely assume it to be Einstein, i.e. Ric = Kg. Its effect is 
essentially an additive constant (related to A) in all equations and constraints. This 
situation does not seem to be systematically investigated in the literature. In partic- 
ular, one can wonder whether the introduction of A is "catastrophic" or in contrary 
produces only a moderate effect. A similar situation is that of the paradigmatic exam- 
ple in holomrophic dynamics, of the quadratic family z h- > z 2 + c. Here the variation 
of the parameter c generates a chaotic dynamics as well as a fractal geometry Q. 



10.8. Wick rotation. Here M = I x M is endowed with the Riemannian metric 
g = +dt 2 + gt- Writing Ric = 0, yields: 
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(30) 



dt Kt 



-2h 



Ric + tr gt (k t )k t - 2k tqt k 







r + \kf 



at 



(tr gt kt) 2 (Constraint) 



For instance, this allows one to construct examples of Riemannian Ricci flat man- 
ifolds, of co-homogeneity 1, i.e. their isometry group has codimension 1 orbits. 

Notice that the (true) Einstein equations (i.e. without symmetries) can not be 
solved in a Riemannian context (they cannot be transformed to a hyperbolic PDE 
system). Maybe, this Bianchi situation can give insights on the reasons behind this 
fact. 

Finally, it does not seem there exists a "true Wick rotation", i.e. some corre- 
spondence between solutions of Bianchi-Einstein equations in the Lorentzian and 
Riemannian cases. (Compare with Q). 



10.9. Orthonormal frames approach vs Metric approach. As a result of a search 
on fundamental references in this area, "dynamical systems and cosmology", one 
can get at least (71 [19) and ll20l which are surely the most known and recent synthesis 
in this "emerging" domain. The authors adopted there an "orthonormal frames ap- 
proach" in opposite to our "metric approach" here (see explanations therein). They 
obtained the following system of quadratic polynomial differential equations on R 5 . 



(31) 



N> 
N 2 



-(2 

-(2-g)S_ -<S_ 
(q - 4S+)iV 1 
(g + 2S+ + 2V3S_)iV 2 
(g + 2S+ -2V3S_)iV 3 



where: 



S- 



(N 2 - N 3 ) 2 - Nt (2Nt -N 2 - N 3 ) 



^= (N 3 - N 2 ) (iVi - N 2 - N 3 ) 



q = i(3 7 -2)(l-K) + |(2- 7 )(SUS 2 _) 

K = ^[Nl + N 2 + N 2 -2{N l N 2 + N 2 N 3 + N 3 N 1 



Here 7 is a parameter: 2/3 < 7 < 2. (See for instance (22)). 
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10.9.1. Comparison. This system of differential equations must be "equivalent" to 
our equations (£ jl0.2l) on the tangent bundle of a Milnor flat, which was a rational 
differential system on M 6 with one constraint. A formal definition of equivalency 
of approaches is that the two systems are "bi-rationally equivalent". However, the 
transformation of our system to this polynomial system is by no means obvious. 
This last system was not a priori motivated by simplifying our more "naive" one, 
but rather by considering another point of view in considering Einstein equations. 
Instead of studying the evolution with time of the metrics on spacelike slices, one 
considers the evolution of brackets of orthonormal frames on these slices. The gauge 
freedom is more subtle in this case, but still this method is very clever, as shown by 
the simplified form of the equations here. In our Bianchi case, i.e. where spacelike 
slices are Lie groups with left invariant metrics, one can very roughly say that the 
bi -rational equivalence comes from the projection map Mil{Q) — > Sym(Q), where 
Mil(Q) is the space of Milnor bases of Q. The next step is to lift the Einstein equation 
(including a gauge choice) to Mil(Q) (more precisely an associated bundle) and to 
take the quotient by the G-action! 

10.10. Further remarks. This beauty of Sym n appeals one to go beyond..., but 
as we said, our contribution here is essentially preliminary and expository. Let us 
mention some facts that were not considered here (with the hope to give details on 
some of them in the future). 

10.10.1. Variants of Sym^. First, one can generalize the discussion from Sym+ 
to Sym* n , the space of all pseudo-Euclidean products, i.e. non-degenerate quadratic 
forms. Everything extends there, a pseudo-Riemannian metric (on Sym* n ), Ricci 
maps, Bianchi-Einstein flows... 

The components of Sym* n are spaces of quadratic forms of a given signature. As 
for Sym+, each component is a pseudo-Riemannian symmetric space and plays a 
universal role in its class. 

• Complex case. The same is true for complex spaces: Sym* n (C), the space of 
complex non-degenerate quadratic forms on C n , is a holomorphic symmetric space... 

•Projectivization. Taking the associated projective spaces will send all these spaces 
into compact ones, and hence compactify them, by attaching various boundaries, with 
more or less nice interpretations. A natural requirement is that ideal points corre- 
spond to collapsing of Riemannian metrics, say in the Gromov sense ifTTTl (restricted 
here to homogeneous spaces). By algebraicity, all differential equations extend to the 
projective spaces. 
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• Fiberwise constructions. If E — > B is a vector fiber bundle, then one can asso- 
ciate to it Sym + (E)... 

• Configuration spaces. Another interesting aspect of Sym+ is its configuration 
space aspect. We mention here the case of " hydrodynamics", where a geometric 
formalism (a Riemannian metric, its geodesic flow...) was developed (see for instance 
0) following similar ideas as those presented here. There are also other non-linear 
and infinite dimensional situations, in particular a space Sym* (E) associated to a 
HUbert space E could be exciting! 

10.10.2. Geodesic flows of left invariant metrics. For a Lie group G, Sym + (Q) 
plays a role of a parameter space of its left invariant metrics. The geodesic flow 
of any such metric on G is a second order quadratic ODE system on Q (H. It is 
intersting to study the dependence on parameter of the qualitative properties of these 
geodesic flows. 

10.10.3. (Locally) Homogeneous, but "non-simply homogeneous" spaces. Instead 
of left invariant metrics on Lie groups, one can consider general homogeneous spaces, 
say, those endowed with an isometric transitive, but not necessarily free action of 
a given group G. More important is the case of locally homogeneous spaces, i.e. 
when the metric varies in the space of all those locally modeled on a fixed space X 
endowed with a (non-fixed) G-invariant metric (but G is fixed). Here, G does not 
act (it acts only locally, as a pseudogroup). As an example, we have the Robertson- 
Walker-Friedman-Lemaitre spacetimes lTT2l[T8ll24l . which are warped products M = 
T x w N, g = —dt 2 + w(t)g, where N has a constant curvature. 

10. 10.4. Dimension 2+1. So far, only the Gauss gauge has been considered. Maybe, 
this is because of its "deterministic character", i.e. it gives rise to autonomous dif- 
ferential equations, instead of non-autonomous ones, as in the generic case. There 
are however other situations where interesting gauges are available. As an example, 
in 't Hooft's theory of systems of particles in dimension 2 + 1 ll23l . one has a flat 
polyhedral surface with singularities, evolving (locally) in a Minkowski space. The 
gauge here is fixed by the fact that time is locally equivalent to a "linear time" in the 
Minkowski space. In particular the time levels remain flat polyhedral. By considera- 
tion of suitable spaces of such surfaces, one may be convinced there is a configuration 
space approach similar to our situation here. 

10.10.5. Non-empty spaces. Recall that M = I x M is a perfect fluid, if Ric = 
(p + p)dt 2 + pg, where p is the pressure and p is the density |[T2llT8ll2"4"l . A Bianchi- 
Einstein flow can be defined in this case, when p and p are functions on TSym+, 
or more reasonably, when they are basic functions, i.e. they depend on the coordi- 
nate q € Sym^ alone. Robertson-Walker spacetimes are examples of perfect fluids 
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(strictly speaking, they would be covered by our approach, once we consider general 
locally homogeneous spaces, as discussed above). 

10.10.6. Quantization of the Bianchi-Einstein flow. We strongly believe this is a nat- 
ural case that can be treated by a quantum gravity theory (see for instance [8]), that 
is, a reasonable quantization of the Bianchi-Einstein flow should be possible... 

10.10.7. A modified Einstein equation. World would be perhaps simpler if the Ein- 
stein equation on TSym^ were given by the mechanical system determined by the 
Riemannian metric (, ) on Sym+ as a kinetic energy, and the Hilbert action % as a 
potential energy. Recall [1] that solutions of such a mechanical system are curves 
q(t) G Sym^, satisfying 

Vg>(t)J(t) = -Vft(<?(t)) 

(V is the Riemannian-connection and VH is the Riemannian gradient of H). Other 
more "realistic" modified equations are obtained by replacing the Riemannian metric 
by a pseudo-Riemannian or a Finsler one of the form L a p or F a $ (©. 
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